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It is argued that many non-extremal black branes exhibit a classical Gregory-Laflamme 
fN| , instability if, and only if, they are locally thermodynamically unstable. For some black 

branes, the Gregory-Laflamme instability must therefore disappear near extremality. For 
the black p-branes of the type II supergravity theories, the Gregory-Laflamme instability 
disappears near extremality for p = 1, 2, 4 but persists all the way down to extremality for 
p = 5,6 (the black D3-brane is not covered by the analysis of this paper). This implies 
■ that the instability also vanishes for the near-extremal black M2 and M5-brane solutions. 

o : 



Abstract 



1 Introduction 

^ 1.1 The Gregory-Laflamme instability 

• i-H ! The Gregory-Laflamme (GL) instability is a classical instability of non-extremal black branes 
[|. The existence of this instability is often explained by arguing that a black p-brane has 
a lower entropy per unit p-volume than an array of black holes with the same total mass and 
charge per unit volume^, and must therefore be unstable. However, this explanation leaves 
something to be desired because it uses thermodynamic (and therefore quantum) arguments to 
explain a classical phenomenon. 

It is, of course, well-known that the classical laws of black hole mechanics || can be re- 
interpreted as thermodynamic laws once it is appreciated that black holes emit thermal ra- 
diation [|J and therefore have a temperature and associated entropy. Conversely, the laws of 
thermodynamics applied to black holes can sometimes be reformulated in classical language. 
The entropy argument above would become an argument based on the second law of black hole 
mechanics which states that the area of the event horizon is non-decreasing in any physical 
process. Since the area of the horizon would increase if a black brane were to evolve to an array 

1 Assuming that the black holes can carry the charge of the black brane. If not, then one might instead 
consider an array of uncharged black holes threaded by an extremal black brane. 
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of black holes, the second law certainly permits a classical instability of black branes. However, 
it does not require it. 



1.2 The Gubser-Mitra conjecture 

The GL instability has recently been re-examined by Gubser and Mitra (GM) || [7|] from the 
perspective of the AdS/CFT correspondence ||. It was argued that the absence of tachyonic 
modes in the AdS glueball calculations || |l(| is evidence that near-extremal D3-branes are in 
fact classically stable, and that this is consistent with the work of GL because GL only checked 
that an instability exists for branes that are far from extremality. However, this appears to 
contradicts the claim of GL JIT) that their instability persists all the way down to extremality^. 
Furthermore, to apply AdS/CFT arguments one must first take the decoupling limit || of the 
black brane, which makes it only infinitesimally non-extrema^. 

More interestingly, the following conjecture was made by GM: a black brane with a non- 
compact translational symmetry is classically stable if, and only if, it is locally thermodynamically 
stable. The definition of local thermodynamic stability adopted by GM was the one used in 
1"2|| , namely that the Hessian of the entropy with respect to the other extensive thermodynamic 



variables (i.e. mass and charges) should be negative definite. This is the mathematical statement 
of the condition that it shouldn't be possible for the system to gain entropy by splitting into 
several parts whilst conserving energy and charge, i.e., it is a statement about the microcanonical 
ensemble. 



At first sight, the GM conjecture appears to contradict the claim [ITJ that the GL instability 
persists all the way down to extremality. For example, near-extremal D3-branes have positive 
specific heat, which implies that they satisfy the above definition of local thermodynamic sta- 
bility and should therefore be classically stable according to the conjecture. However, the work 
of GL only covered black branes carrying magnetic charge with respect to a Neveu-Schwarz 
n-form. Therefore, the only black brane of type II supergravity covered by their work is the 
NS5-brane. This brane does indeed have negative specific heat all the way down to extremality, 
so the fact that it is classically unstable agrees with the GM conjecture. One might think that 
the work of GL could be extended to the F-string by dualizing the NS 3-form, but this changes 
the coupling to the dilaton, taking one outside the class of theories considered by GL. One of 
the results of the present paper will be to show that the GL instability occurs for many more 
solutions than just the ones considered by GL. 

The condition that there should be some non-compact translational symmetry was included 
in the GM conjecture to exclude black holes from the scope of the conjecture. This is because 
it is well-known that many black holes, for example Schwarzschild, are classically stable |H| yet 
have negative specific heat and are therefore locally thermodynamically unstable. 

The GM conjecture can be regarded as a refinement of the entropy argument mentioned 
in the first paragraph of this paper. This entropy argument refers to global thermodynamic 
stability. It seems much more plausible that, if there is a connection between classical and 



2 That is, for any non-zero value of the non-extremality parameter. It was argued that the instability is 
absent for extremal branes, which is as one would expect in a supersymmetric theory. 

3 Thanks to Simon Ross for emphasizing this point to me. Throughout this paper, "near-extremal" will mean 
that some non-extremality parameter lies in a finite neighbourhood of its extremal value. 
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thermodynamic stability, it should be local thermodynamic stability that is relevant. For ex- 
ample, one might imagine some kind of potential barrier separating the black brane from an 
array of black holes, so the only way that the brane could gain in entropy would be by quantum 
tunneling through this barrier. In that case, since the instability would be quantum mechanical 
in nature, global thermodynamic instability would not imply classical instability. 

The evidence that GM gave for their conjecture was based on studying black holes solu- 
tions of gauged M = 8 supergravity in four dimensions. These solutions are asymptotic 



to anti-de Sitter space and the black holes considered were large compared with the AdS ra- 
dius of curvature. This presumably implies that the curvature of their horizons should not be 
significant so they would be expected to exhibit similar behaviour to black branes. Therefore 
if the conjecture is correct, then it should apply to these black holes too. It was found that 
a classical instability appeared almost precisely where local thermodynamic stability was lost, 
with a small discrepancy arising from numerical error. Note that large asymptotically flat black 
holes would not be expected to behave like black branes because there is no length scale with 
which the curvature of the horizon can be compared. 

1.3 Proving the conjecture 

In this paper, a semi-classical proof of the GM conjecture will be given, using the Euclidean 
path integral approach to quantum gravity. This has been succesful in providing a semi-classical 
explanation of other connections between classical and thermodynamic properties of black holes, 
for example the relationship between the entropy of a black hole and the area of its horizon 
15fl . It is convenient to review some basic facts about this approach |L6]. The Euclidean path 



integral for the canonical ensemble is 

'd[g]exp(-/[g]), (1.1) 



where I denotes the Euclidean Einstein-Hilbert action. The path integral is taken over Rie- 
mannian manifolds (M, g) that are asymptotically flat, i.e., tend to the flat metric on S 1 x R 3 
outside some compact region. The S 1 should have proper length (3 = 1/T, where T is the 
temperature. 

This path-integral is only well-defined in the semi-classical approximation. In this approx- 
imation, one looks for saddle-points of the Euclidean action, i.e., solutions of the Euclidean 
Einstein equations. These are referred to as gravitational instantons. Near a gravitational 
instanton, the metric can be written as 

9^u = g^v + Sg^ u , (1.2) 

where g^ u is the metric of the instanton, and 5g^ v denotes the perturbation. One can add gauge 
fixing and ghost terms to the action to deal with the gauge freedom in 5g^ u . The action can be 
expanded around the saddle point: 

/[g] = / [g]+/ 2 [g,5g], (1.3) 
where I2 is quadratic in the fluctuation. The contribution to I2 from the trace of the metric 



perturbation is a wrong-sign kinetic term. This is the well-known conformal factor problem [17| . 
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Nothing new will be said about this problem in the present paper. The part of I2 involving the 
traceless part of the metric perturbation is, in transverse gauge, proportional to 

J d A x^H^A L H^ (1.4) 

where is transverse traceless and A^ is the Euclidean Lichnerowicz operator. The path- 
integral over the transverse traceless part of the metric perturbation can be evaluated by ex- 
panding in eigenf unctions of Al, i.e., solutions of 

A L H„ V = XHfj, u . (1.5) 

The boundary conditions are that should vanish fast enough at infinity to have finite norm, 
and be regular everywhere else. The contribution to the partition function from the instanton 
is 

Aexp(-J [g]) ) (1.6) 

where A is a functional determinant that arises from the path-integrals over the metric per- 
turbation and ghosts. In particular, it includes a factor of (det Al)~ 1//2 . This implies that if 
Al has a negative eigenvalue then there will be an imaginary contribution to the path integral. 
This pathology arises from the fact that the canonical ensemble for gravity is not well-defined 



in infinite volume [18, 16 1, i.e., it is an indication of a thermodynamic instability. 



For the Euclidean Schwarzschild solution, there is a single transverse traceless negative mode 



19|, 20], i.e., a solution to equation 1.5 with negative A. The eigenvalue is [20 



A = -0.19M" 2 , (1.7) 

where M denotes the Schwarzschild mass. The existence of this negative mode implies that 
the Euclidean Schwarzschild solution is merely a saddle point of the Euclidean action, rather 



than a local minimum [19|. Therefore this instanton cannot be regarded as giving an approx- 



imation to the partition function. The correct interpretation is as an instanton describing a 



non-perturbative instability of flat space at finite temperature |20 |. 

It is now possible to explain the main idea of this paper. Consider an uncharged black 
p-brane in 4 + p dimensions with metric 

ds 2 = g liv dx^dx v + 5 ij dz i dz j , (1.8) 

where denotes the four dimensional Lorentzian Schwarzschild metric and z % denote the flat 
spatial worldvolume directions of the brane. GL considered |I| metric perturbations of this 
solution of the form 

h^ v = exp (ifiiZ l )H^ u , = h {j = 0, (1.9) 

where transverse traceless gauge is used for h^, and H^ u is spherically symmetric. GL looked 
for modes that are regular on the future horizon of the brane and decay at infinity, but grow 
exponentially in time: H^ u oc exp (fit). The equation of motion for the perturbation is [[J 

A L H^ = -fM 2 H^, (1.10) 

where denotes the four dimensional Lorentzian Lichnerowicz operator and p? = Y^i Hif^i- 
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GL showed that there is a critical value /x* such that for every //<//*, equation |1 . 10 



has a solution that is regular on the future horizon and decays at spatial infinity, but grows 
exponentially in time. For \i > /x*, no growing mode was found and the solutions are oscillatory 
in time. It follows that the solution with /x = /x* must be independent of time. This solution 
will be referred to as the threshold unstable mode. Note that the classical second law implies 
that a black brane compactified on a sufficiently small torus cannot evolve to an array of black 
holes, which suggests that the black brane should be stable against small wavelength (large /x) 
perturbations, and therefore that a time-independent mode must exist if the brane is unstable 
to large wavelength perturbations. So one can argue that the existence of such a mode is 
necessary for a classical instability without referring to the numerical results of [|l|, ^ . This will 
be important for the more general black brane solutions considered later in this paper. 

It is now possible to make contact with the Euclidean path integral by Wick rotating H^ u 
for the threshold unstable mode. Since is static and spherically symmetric, Wick rotation 
simply changes the sign of H tt - H^ v is translationally invariant on the brane world- volume so 
it can be dimensionally reduced to four dimensions. Wick rotation transforms A^ into A^, 
so equation |1.10| turns into equation [L5] with A = — /x^. Hence it appears that the classical 
Lorentzian threshold unstable mode corresponds precisely to a Euclidean negative mode. It will 
be shown later that the Wick rotated mode does indeed satisfy the correct boundary conditions 
for a negative mode. 

This argument can be checked by comparing the numerical results of |^0| and [|TJ . The latter 
gave results for horizon radius r + = 2 [0, which corresponds to M = 1, i.e., A = —0.19. This 
would give /x* = 0.44, which is in excellent agreement with figure 1 of ffl. It is possible to extend 



this matching to the other dimensions for which the negative mode has been calculated |2T 
namely A = — 1.60r^ 2 in five dimensions and A = — 3.39r^ 2 in seven dimensions. For r + = 2, 
this gives /x* = 0.63 and /x* = 0.92 respectively, again in good agreement with [[Lp. A simple 
scaling argument shows that this agreement must extend to all values of r + . 

It is clear that the classical GL instability of the black brane is very closely related to 
the thermodynamic instability of the Schwarzschild black hole. The threshold unstable mode 
is related to the Schwarzschild negative mode by discarding the exp(i/Xj;z*) factor and Wick 
rotating^. The existence of negative modes is closely linked to local thermodynamic stability, as 
will be explained shortly. It therefore seems likely that the above argument can be generalized 
into a proof of the GM conjecture for more general black brane solutions. That will be the 
subject of the rest of this paper. 



1.4 Outline of this paper 

The next section describes the class of black branes that will be discussed in this paper, namely 
magnetically charged dilatonic black branes. In section 0, the relationship between local thermo- 
dynamic stability and negative modes is explained. Section | discusses classical perturbations 
and in section [5] it is shown how these are related to negative modes. The final section summa- 
rizes the main conclusions and suggests directions for future research. 

4 This argument shows how a negative mode of the Schwarzschild black hole is related to a static zero mode 
of the black string. This relationship was noted in jjffi] , although the connection with classical stability was not 
appreciated (indeed it was argued that the black string is classically stable!) 
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2 Black brane solutions 



This paper will discuss black p-brane solutions of the action 

I = eJ d D x^g (-R + 1(«90) 2 + ^Ff n 



) ' 



(2.1) 



where e = +1 for Euclidean signature and —1 for Lorentzian signature. Fi n ) denotes a n-form 
field strength obtained from a (n — l)-form potential A( n _i) by F( n ) = cL4( n _i). 



Black brane solutions of this action have been given by several authors 123L 24, 25, 26 



The notation of [26| will be used here. The magnetically charged black brane solutions can be 
written: 



e 



1 H = sinh 2 a 



+ 1 H = sinh u 

^ i 2 

= 1 H — j sinh /x, 



4d 

' A(D-2) 



A(D-2) 



ee 2/ dt 2 + dz l dz' 



(e- 2f dr 2 + r 2 d{l 



2 



,2/ 



1 ^ 



(2.2) 



The parameters d, d and A are defined by 

d = n — 1, d = p + l = D — d — 2, 



A = a 2 + 



2dd 

d-: 



(2.3) 



The n-form field strength is a constant times the volume form on S n . The two free parameters 
k and \x are related to the charge A and mass per unit p-volume m by p6| 



A 



dk 



The temperature is 



sinh2/i, 



T 



m = k I — sinh /x + a + 1 



47rr i 



(cosh /x) 4//A 



(2.4) 



(2.5) 



The extremal limit is k — > 0, /x — > oo with /ce 2M fixed. It is straightforward to verify that the 
specific heat at fixed charge is always negative if 2d/A < 1. However, if 2d/A > 1 then this 
specific heat is negative far from extremality but positive near extremality. 

For the black p-branes of the type II supergravity theories J23], with the exception of the 
D3-brane, one can always dualize to make the solution magnetically charged. The solutions 
then take the above form, and one finds that the specific heat changes sign near extremality 
when p < 4 (i.e., DO, Fl, Dl, D2, D4) but not for p = 5, 6 (NS5, D5, D6) (Dp-branes with 
p > 6 will not be considered here because they are not asymptotically flat). For the solutions 
considered by Gregory and Laflamme, the specific heat never changes sign. The specific heat 
of the D3-brane does change sign but this brane lies outside the class of solutions considered 
in this paper. Nevertheless, it seems likely that the arguments of the present paper could 
be generalized to include this case. For the black M2 and M5 brane solutions |24| of eleven 
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dimensional supergravity, the specific heat is negative far from extremality but becomes positive 
near extremality. 

Dimensional reduction on the worldvolume of the above solution leads to a magnetically 
charged dilatonic black hole solution in n + 2 dimensions. The properties of such solutions were 



extensively discussed in [27]. It is clear that the temperature of this black hole must be the 



same as the temperature of the black brane. It is also the case that the charge of the black hole 
is the same as the charge of the black brane, and the mass of the black hole is proportional to 



the mass per unit p- volume of the black brane [281. It follows that the specific heat at constant 



charge of the black hole has the same sign as the specific heat at constant charge of the black 
brane. 



3 Negative modes and the specific heat 
3.1 Canonical ensemble 

The canonical ensemble for gravity can be rendered well-defined by imposing boundary condi- 
tions at a finite boundary rather than at infinity |29|]. Consider a spherical cavity of surface 



area 47rr^, whose boundary is held at a fixed temperature T. The partition function in the 
canonical ensemble is given by a path integral over Riemannian geometries (M, g) such that 
DM = S 1 x S 2 , where the S l has proper length /3 = 1/T and the S 2 has area Attt"^. 

In the case of a pure gravity theory, for small Trg the only solution to the Euclidean Einstein 
equations that obeys these boundary conditions is flat space with a periodically identified time 
coordinate. For larger Trg, there are three solutions: flat space, a small Schwarzschild black 
hole, and a large Schwarzschild black hole (with size of order r#). The smaller black hole has 
negative specific heat and the larger black hole has positive specific heat. 

It is straightforward to demonstrate that the black hole with negative specific heat cannot 
be a local minimum of the Euclidean action and must therefore have a negative mode. For a 
Schwarzschild black hole, this was shown in [001 using a construction given in |31[j . The idea is 
to construct a 1 parameter family of geometries for which the Euclidean action (relative to flat 
space) takes the form 

J = 0E(x) - S(x), (3.1) 

where x is the parameter labelling the path. There should be should be some value x = x(T) 
for which the geometry becomes the black hole geometry of interest. At this value, E and S 
are the energy and entropy of the black hole: 

E(x(T)) = E(T), S(x(T)) = S(T). (3.2) 

For other values of x, the interpretation of the functions E and S is not important. The black 
hole extremizes the action, so at x = x(T), 

which is simply the condition for thermal equilibrium of the black hole with a heat bath of 
temperature T: 



AE dS 



dE 



The second derivative of the action along the path parametrized by x is 



d 2 I\ („dx\ 2 dE 




dx 2 ) T \ dT) dT 



(3.5) 



It is evident from this equation that if the specific heat dE/dT is negative then the black hole 
cannot be a local minimum of the action and therefore it must have a negative modeQ 

In order to exhibit this connection between specific heat and negative modes it was necessary 
to assume that there exists a 1-parameter family of geometries, including the black hole, for 
which equation |3J] is obeyed. This will now be demonstrated for the black holes of relevance 
to this paper, namely those obtained by dimensional reduction of the black brane solutions 
given above. These are dilatonic black holes in n + 2 dimensions carrying magnetic charge with 



respect to a n-form field strength. The construction is a straightforward generalization of pi 
The Euclidean metric can be written 

ds 2 = U{r)dr 2 + Vir^dr 2 + R{r) 2 dQ 2 n . (3.6) 

The exact form of U(r), V(r) and R(r) will not be important. The dilaton is a function of r 
alone, and the n-form is proportional to the volume form on S n . 

In the canonical ensemble, if the boundary is at r = then the proper length of the 
Euclidean time direction at r = r# should be (3 — 1/T. It is convenient to normalize r so that 
t ~ t + /3, which implies U{tb) = 1- The functions U(r) and V(r) both vanish at the horizon, 
r = r + . The condition for regularity at r = r + is 

47T 



sJU'(r + )V'(r + ) = j, (3.7) 

which gives r + = r + (T). The black hole solution is specified by the parameters k and /i, which 
are determined by the temperature T and charge A. This black hole will be referred to as 
B(T,X). 

The other geometries in the family will be taken to have a metric of the form |3.6| but now 



with U(r) treated as an arbitrary function. V(r), R(r), the dilaton and the n-form are all taken 
to be the same as for a black hole with charge A but arbitrary temperature. In other words, 
k and \i are constrained only to keep A fixed but not T, so this gives a 1-parameter family 
of metrics on the surfaces of constant r. This family can be parametrized by the value of r + . 
Choosing r + = r+(T) reproduces the geometry of the spatial sections of B(T, A). 

This choice of geometry for surfaces of constant r implies that the constraint equations of 
general relativity are satisfied on these surfaces, r is still taken to be identified with period f3. 
The only restrictions on U(r) are U(tb) = 1 (so the metric obeys the boundary conditions), 
U(r + ) = (so the topology is the same as that of the black hole B(T,X)), that U'(r+) obeys 
equation |3.7| (so that the metric is regular at r = r+) and that U(r) is positive for r > r + . It 
is obviously possible to choose a suitable family of functions U(r) so that the metric becomes 
that of B(T, A) when r + = r + {T). 



5 Note that the path of off-shell geometries parametrized by x does not correspond directly to an eigenfunction 
of the Lichncrowicz operator, but rather to some linear combination of eigenfunctions. Since the action decreases 
along this path, at least one of these eigenfunctions must have a negative eigenvalue. 
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The action of these geometries (relative to flat space) can be calculated using the methods 



of |I6|, f32j. One first excises the region r < r + + e. In the region r > r + + e, one write the action 



in Hamiltonian form. The contribution from this region is then simply 

1+ = (3H, (3.8) 

where H is the Hamiltonian evaluated on the above metric. Since the above metric is static 
and satisfies the constraint equations, the only contribution to H comes from a surface term at 
r = tb-i giving 

I+ = PE, (3.9) 

where E depends only on r + and the charge A. The total action is given by adding the action 
of the region r < r + + e. In the limit e — > 0, this is given entirely by a gravitational surface 
termf] at r = r + which, after using equation |3]7] reduces to 



A 

1g' 



(3.10) 



where A is the area of the horizon at r = r + , and G is Newton's constant. The total action is 
therefore 

I = I++I- = PE-S, (3.11) 

where S = A /AG, which depends only on r + . This establishes equation |3.1| (with x = r + ). 
The proof that negative specific heat implies the existence of a negative mode then proceeds 
as above. The only subtlety is that, since the charge was held fixed above, the specific heat 
should be calculated at constant charge. The reason for working in the canonical ensemble 
(fixed charge) rather than the grand canonical ensemble (fixed potential) is that the theory 
does not contain charged particles and therefore the charge cannot fluctuate. This is also true 
for electrically charged black holesQ 



3.2 Physical interpretation 

There is a simple physical interpretation of the off-shell geometries considered above. Consider a 
typical member of this family. The geometry of the spatial sections is simply that of a black hole 
of charge A but with temperature T' ^ T. Since the metric is static (independent of r) analytic 
continuation to Lorentzian signature at a surface of constant r yields initial data for solving 
the Lorentzian Einstein equations and the unique solution is the black hole with charge A and 
temperature T' . In other words, the path of off-shell geometries constructed above corresponds 
to an instability in which the black hole changes temperature at fixed charge, which is precisely 
the instability that one would expect for a black hole in contact with a heat bath. Note that the 
off-shell geometry is not simply given by analytic continuation of the metric of a black hole of 
temperature T' since this would not give a metric with the correct period for the Euclidean time 
at r = rs- Note also, that if this interpretation is correct then many different U(r) correspond 
to the same Lorentzian black hole. 

6 A surface term involving the n-form vanishes for magnetic black holes. 

7 In the case of electric charge (for n = 2), it is necessary to add an extra surface term to the Euclidean 
action when working in the canonical ensemble fl33L pij . 
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It has been shown that negative specific heat implies the existence of a negative mode. 
To see that this is a physical negative mode (rather than a negative mode arising from the 
conformal factor problem), consider the off-shell geometries near B(T, A). These can be regarded 
as perturbations of B(T, A). The function U(r) is essentially arbitrary for r + < r < tb, so by 
adjusting U(r), the trace of the metric perturbation can be freely adjusted. If the negative 
mode were associated with the conformal factor problem then one would expect increasing 
the trace of the metric perturbation to make the action more negative. However, since the 
constraint equations are satisfied on surfaces of constant r, the action does not depend on U(r) 
and therefore this change in the trace of the metric perturbation does not change the action. 
It follows that the negative mode whose existence was demonstrated above is a non-conformal 
negative mode. 

The above discussion considered black holes in finite cavities. This is a rather unphysical 
situation since it is hard to see how an appropriate cavity could be constructed. The rest of this 
paper will consider only black holes in infinite space, i.e., the tb — > oo limit will be assumed. 

3.3 Other ensembles 

Since the charge of the black hole cannot fluctuate, the condition for local thermodynamic 
stability in the microcanonical ensemble is that the second derivative of the entropy with respect 
to the mass must be negative. Equivalently, the specific heat should be positive. So the 
condition for local thermodynamic stability of the black hole is the same in the canonical and 
microcanonical ensembles. 

In a theory containing charged particles, one might choose to use the grand ensemble for 
electrically charged holes. This could lead to new negative modes involving the gauge field. 
For example, consider the Reissner-Nordstrom solution. The necessary and sufficient conditions 
for local thermodynamic stability in the grand ensemble is that the Hessian of the Gibbs free 
energy with respect to the entropy S and charge Q be positive definite. This is equivalent to^ 

(dE\ ( d<&\ 

where $ denotes the potential at infinity (or at r = rs) in a gauge regular at the horizon. 
In four dimensions, one finds that the former equation is violated for Q/M < v3/2 but the 
latter equation is violated for Q/M > Vo/2. Thus one would expect an instability analagous to 
the thermodynamic instability of the Schwarzschild black hole for small Q, and the associated 
negative mode would involve the metric as above. However, for large Q, the instability involves 
the black hole gaining or losing charge, so the instability would involve the electromagnetic 
field and charged matter fields, as for the black holes studied in || [7|. Hence, in this case one 
would expect a negative mode involving the electromagnetic and matter fields. In a theory that 
did not contain charged matter, one would work in the canonical ensemble and the second of 
the above stability conditions would not be required so one would expect no negative mode for 
Q/M > a/3/2. This is confirmed by numerical calculations which show that the negative 
mode involving the metric vanishes precisely at Q/M = Vo/2. 

8 See section 21 of [B5j and replace — V with Q and p with 
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Note that it has only been shown that negative specific heat implies the existence of a 
negative mode. Conversely, the existence of a negative mode implies that the black hole is only 
a saddle point of the Euclidean action, which implies that it suffers from some pathology in the 
canonical ensemble. This suggests that it should be locally thermodynamically unstable. The 
Reissner-Nordstrom example is evidence that this is true, but I am unaware of a more general 
proof. 



4 Classical stability 

The aim of this section is to derive the equations of motion for classical perturbations of the 
black brane solutions discussed in section |[ 

4.1 Equations of motion 

The analysis will be simplified by a convenient choice of conformal frame. For the moment the 
choice of conformal frame will be kept arbitrary by writing 

9mn = e 2b *g M N, (4.1) 
where a caret refers to the D-dimensional Einstein metric. In the new frame, the action is 



I = ef d D x^g [e-» (-R + k{d<P) 2 ) + ^e^if n) 



(4.2) 



where 



/? = -(£> -2)6, k = i - (D- l)(D-2)b 2 , a = a + (D - 2n)b. (4.3) 
The equations of motion can be written 

V M (e^F MN '- N ^) = 0, (4.4) 

V 2 - /?(<90) 2 = ^ a+ ^F\ (4.5) 
2n\ 

Rmn = (k + (3 2 )d M d N( j)-f3VMVN(J> 

+ ^^^Fm^^F^^ - A e (« + W FW (4.6) 

1 D-2 J \ D-2 } D — 2 K ' 



where 



The conformal frame that will be used is given by 



d 

a(D-2)' 
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In this frame, there is no conformal factor multiplying the flat worldvolume directions of the 
black brane metric. Making this choice has the added advantage that it gives 



A = 0. (4.9) 

Perturbing the equations of motion gives 

V N (e«*SF NPl ~ p »- 1 ) + {ad N 8(f>)e a ^F NPl - p ^ - {n - l)e a *F NQ ^- p ^V N h Pl] 

~h N M V N (e^pMPx.-P^ _ Vn ( h N _ KgN\ e afpMP 1 ...P n . 1 = q (4 1Q) 



V 2 8<b - h MN V M V N 4> + ph MN d M (f>d N (f) - d M ^V N (h MN - \hg MN ^j - 2(380 ■ 



—> +f3) \a + (3)5(bF 2 = 0. (4.11) 
2n\ 



V 2 h M N - Vm V p ( h PN - -hg PN ) - V N V P (h PM - -hg 



2 a J V 2 

1 



PM 



2-Rp(M^7V) — 2R(M\PQ\N)h PQ + 2(3 ^V(M^) — -V p h,MN ) 9f 

2(3V M V N 54> + A(k + (3 2 )d {M 0d N 



+ 1 c (*+m 



n 1)t ^(M^^V,..^, - (n - l)h p< ^F MPRl ... Rn _ 2 F NQ R - R — 

+ (a + /3)50F A/Pl ... P „_ 1 F w Pl - p "- 1 ] = 0, (4.12) 

where bars denote background quantities, and the background metric is used to define covariant 
derivatives and raise and lower indices. 

It is convenient to split the D-dimensional indices M,N into indices i,j along the spatial 
worldvolume directions and indices fi, v in the other directions. Indices m, n will be used to 
denote directions on the n-sphere. The choice of conformal frame made above implies that 
connection components and curvature components vanish if they have an % index, z will be used 
to denote spatial worldvolume coordinates and x to denote other coordinates. 

4.2 The Gregory-Laflamme ansatz 

Following Gregory and Laflamme ||, only s-wave perturbations will be considered. For such 
perturbations, 

h tm = h rm = h im = 0, h™ = K(t,r)5™. (4.13) 
Other components, and the dilaton perturbation, are also taken to be independent of the co- 



ordinates on the sphere. This implies that the terms involving Hmn arid 8(b in equation [4.10 
vanish, just as in 0. It is therefore consistent to set SF = 0. In [[|, it was shown that there 
was no unstable mode associated with SF. It seems very likely that the same is true here, so 
5F = will be assumed from now on. 
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GL made the ansatz 

h MN = exp(ifiiZ l )H M N(x), 5<p = exp(ifiiZ l )f(x) (4.14) 
and considered only longitudinal perturbations, i.e., 

H ilt = mA^x), (4.15) 

Hij = HiiijH{x). (4.16) 

The same ansatz will be adopted here. Under gauge transformations, the metric perturbation 
transforms as 

hMN h' MN = h MN — Vm£jv — Vat^a/, (4-17) 
so under a gauge transformation with £m = exp(ijt/j2*)Hjyj and Sj = ifaE the metric perturbation 
becomes 

A; = A.-d.E-E^, (4.18) 
H' = H + 2E. 

The quantity H pv — 2V ^A^ — V M V„if is gauge invariant under such transformations. Note 
that it is possible to transform away A^ and if by a gauge transformation. The following gauge 
choice can therefore be made 

= = 0. (4.19) 

This is not the same gauge as used by Gregory and Laflamme. However, in many respects, 
it is much simpler to use this gauge that the de Donder gauge that they adopted. The main 
reason for using this gauge is to relate the classical instability to a negative mode. The negative 
mode can be viewed as arising from a lower dimensional black hole solution, and is therefore 
independent of the coordinates z\ It therefore appears simplest to adopt a gauge in which the 
perturbation is independent of these coordinates. 

Using the gauge |4.19| and the above ansatz for the form of the perturbation, the dilaton 
perturbation equation reduces to 

-V 2 / + 2W-5/ + FV^ 



- (3H^d^dA + d^V u H» v - -mg 



-hi/ 



2" p " 

- 2{n <X _ 1) f {a+ ^ {H^F, pi ... Pn _ 1 F u ^--^ - ^F 2 f^j = A/, (4.20) 

where^ A = — fi 2 . The reason for writing the equations in this form will soon become apparent. 
The \xv components of the metric perturbation equation are 

- V 2 H^ + 2V {i y p H v)p - V M V„#; + 2R p(ll H p u) + 2R^ pa \ u) H p ° 

- (3 (2V {P H P U) - V p H^) + 2(3V^J - A(k + (3 2 )d { ^d v) f 

+ J^iy e(a+ ^ ( (n " ^ W.a„- 2 ^ Ai --- A "- 2 - F^^F^-^ia + /?)/) 

= XH pu . (4.21) 



9 Note that /i 2 — X^MiMii which should not be confused with the parameter /i that appears in the black 
brane solution. 
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The \ii and ij components of the metric perturbation equation reduce to 



y, = v v m - m v M - 2 f( k + P 2 W = o, 



(4.22) 



Z = H p p -2(3f = 0. (4.23) 

It is straightforward to check that the trace of equation |4.21| together with equations [4.22| and 
4.23| implies equation 4.20 . Equation 4.23 can be used to eliminate / from equations 4.21 and 
4.22| and then the problem reduces to solving equations |4.21| and |4.22| . Note that equation [4.22 
looks rather like a gauge condition. However, this cannot be the case because a gauge has 
already been chosen, so equation [4.22| should be regarded as an equation of motion. 

GL looked for solutions with time dependence e Qt . That could be done here too, with 
equation |4.22| used to express two of the four remaining unknown functions (H u , H tr , H rr 
and K) in terms of the other two. These could then be substituted into equation [4.21| to give 
differential equations for these two remaining unknowns. This appears to give a set of equations 
in far fewer unknowns that those studied in |], |2|. This difference arises because the gauge 
choice adopted above is not the same as the one used by GL. 

This route will not be pursued here since the aim of this paper is to exhibit a connection 
between classical instability and Euclidean negative modes. It was argued above that the 
necessary and sufficient condition for a classical GL instability to exist is that there should exist 
a static threshold unstable mode. For such a mode, 



n = H* 



0. 



(4.24) 



4.3 Boundary conditions 

For a threshold unstable mode, the boundary conditions at the horizon of the black brane are 



particularly simple. If the background metric is written in the form |3]6] then the perturbed 
Lorentzian metric can be written 



ds 2 



-U(r)(l + (f>(r, z))dt 2 + V{r)- l {\ + tp(r, z))dr 2 + R(r) 2 (l + k(r, z))dn 2 + dz l dz\ (4.25) 



where <p, ip and k are infinitesimal perturbations with z-dependence exp(i/ijZ 4 ). To examine reg- 
ularity of this perturbation on the future horizon, one can transform to Eddington-Finkelstein |J 
coordinates. One finds that the metric is regular at the horizon if, and only if, the perturbation 
is bounded as r — ► r + and 

<!>(r + ,z)=il>(r + ,z). (4.26) 
At large r, the boundary condition is that the perturbation should decay. 



5 Thermodynamic instability 

In this section, it will be shown that a classical Lorentzian threshold unstable mode exists if, 
and only if, the black hole obtained by dimensional reduction of the black brane has a Euclidean 
negative mode. 
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5.1 Wick rotation 



If a threshold unstable mode exists then it can be Wick rotated to Euclidean signature (this will 
change the sign of h tt ). This gives a solution to the equations describing perturbations around 
the Euclidean signature black brane solution. As discussed in the introduction for the case of 
the Schwarzschild black hole, the perturbation that will be related to the Euclidean negative 
mode will not be the classical perturbation given by Iimn and 5<f) (which are proportional 
to exp(ifiiZ 1 )), but rather the (off-shell) perturbation given by the translationally invariant 
functions Hmn and /. This perturbation satisfies equations [4.20| and [4.2 1| , which look rather 
like eigenvalue equations for small quantum fluctuations around a Euclidean black hole solution. 
Moreover, since A = —p 2 is negative, it appears that this particular fluctuation has a negative 
eigenvalue and is therefore a Euclidean negative mode of the black hole. In this section it will 
be shown that this is indeed the case. This Euclidean perturbation will be referred to as a 
candidate negative mode. 

Since the background solution and the perturbation in question are both translationally 
invariant, it is often convenient not to distinguish between perturbations of the black brane and 
those of the black hole obtained by dimensional reduction. A negative mode of the black hole 
can be regarded as a translationally invariant negative mode of the black brane. In order to 
avoid having to work in two different conformal frames, the black hole will be studied in the 
frame in which its metric is given by simply neglecting the flat directions of the black brane 
metric. Of course, physical quantities such as masses and charges should be calculated in the 
lower dimensional Einstein frame. However, if a negative mode exists in one conformal frame 
then one must exist in all other conformal frames since a change of variable in the Euclidean 
path integral cannot change the number of negative modes. 

It is necessary to show that the candidate negative mode is really a physical negative mode. 
The first thing to check is that it satisfies the appropriate boundary conditions for a quantum 
fluctuation around the black hole. These are that the fluctuation should decay sufficiently 
rapidly, and preserve regularity at r = r + . The first of these follows from the large r behaviour 
of the Lorentzian solution, which is clearly not changed by Wick rotation. The second can be 
seen by Euclideanizing |4.25| (i.e. setting t = —it), dropping the exp(ipiZ l ) factors, and then 



defining a new radial coordinate p near r = r + by 



P=^=^(l + U(r + )). (.1.1) 



y/V'(r + ) 



2 



In the new radial coordinate, the rp part of the metric becomes, near r = r+, 



dp 2 + P 2 ^^^ (1 - Hr + ) + 0(r + )) dr\ (5.2) 



Regularity of the background solutions implies that r ~ r + (3 and equation [O] is satisfied. The 
perturbation only preserves this regularity if it is bounded and satisfies 

V>(r+) = 0(r + ), (5.3) 

and this is clearly satisfied by the Wick rotated mode if it obeys the Lorentzian boundary 
condition |4.26| . 
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5.2 Negative modes 

The candidate negative mode can be written 



89. 



hi/ 



H, 



/IV 1 



$9ui = $9i 



0. 



(5.4) 



H^ v and / satisfy the Wick rotated versions of equations [4.20| and [4.2 1| , which can be written 

as 

G{x) f d D x'0{x,x')3L(x') = ^/ge~^XX(x) (5.5) 



where 



and 



8 2 I 



s 2 i 



X 



09mn 



Q(x,x') 



5<p(x)S<j>(x') 
S 2 I 



&4>{x)8gp Q (x') 
S 2 I 



G 



2 P n 
D^dPQ 



&9mn (x)8<f>(x') 8g MN (x)Sg PQ (x') 



2/1 2/3 2 



2 / 2/3 2 \ 

^9m(p9q)n — jj^2 V "d^ J 9mn9pq 



(5.6) 



(5.7) 



The matrix G = Gjj can be regarded as the metric on the space of fluctuations around a 
solution of the field equations, i.e., if the fluctuation (SgMN, 84>) is denoted schematically as ipi 
then the kinetic term for the fluctuation is G IJ dipi ■ dipj- The presence of G in equation |5.5| is 
required to balance up and down /, J indices. 



Equation |5.5| is the equation governing quadratic fluctuations around the Euclideanized black 
brane solution. So the candidate negative mode satisfies the correct eigenvalue equation and 
the correct boundary conditions for quantum fluctuations around the black brane, or black hole 
(by dimensional reduction). The remaining thing to check is that it really is a negative mode, 
i.e., a fluctuation that decreases the Euclidean action. A mode with eigenvalue A has action 



51 = X I d D x^e-^X. 2 



(5.8) 



where X 2 = X T G X X = G 13 XjXj . Expanding in terms of / and Ha V gives 



X 2 



B 2 n\ 1 

j^)f + PHPf + -H^H l 



2 



2 1 x uMN fr 



P 

D — 2 
AD 



(5.9) 



where Hmn is defined to be the D-dimensional Einstein frame metric perturbation, with Hmn 
its traceless part. It is clear that the metric G is not positive definite, but has a negative 
eigenvalue associated with the trace of the Einstein frame metric perturbation. In other words, 
the kinetic term for this trace has the wrong sign. This is, of course, the conformal factor 
problem [[17]]. 

Equation implies that modes for which X 2 < everywhere are negative modes (i.e. 
decrease the action) if A > 0. These will be regarded as unphysical because they arise from the 
conformal factor problem. The modes with X 2 > everywhere are associated with directions in 
the path integral along which the kinetic terms remain bounded, and therefore if such a mode 
has A < then it represents a physical non-conformal negative mode. 
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The candidate negative mode was obtained from a Lorentzian solution with Z = 0, so it has 



2/3/. This implies 



X 2 



2{3 2 



D 



2JP_\ 2 \ E ^ E , 
n + 2 J 2 ^' 



(5.10) 



where is the traceless part of H^. A sufficient condition for X 2 > is that the coefficient 
of f 2 should be non-negative. Using equations fO| , fO| and [2~5| , it is easy to show that this 
occurs if, and only if, 

, 2(D - p - 3) 2 (2 - p) 

a > —, w — r^- (5.11) 

" (£>-2)(£>-p) 1 J 

This is automatically satisfied for p > 2. For p = 1, the cases of interest have to be considered 
individually. For a Dp-brane, D — 10 and a = (p — 3)/2 and it follows that equation [5.11 
is satisfied for p = 1. For the magnetic NS branes considered in 0, one has -D = 10 and 
a = (7 — p)/2, and equation |5.11| is satisfied for p = 1. Thus, for all the cases of interest, the 
candidate negative mode is a non-conformal negative mode. 



5.3 Proving the converse 



It has been shown that the classical Lorentzian threshold unstable mode can be converted into a 
non-conformal Euclidean negative mode by dropping the exp(i/ijZ*) factors and Wick rotating. 
If one accepts that the existence of a non-conformal negative mode implies local thermodynamic 
instability then it follows that the existence of a Gregory-Laflamme instability implies that there 
is also a local thermodynamic instability. 

To establish equivalence of classical and local thermodynamic stability, it is necessary to 
prove the converse. As shown above, local thermodynamic instability (negative specific heat) 
implies the existence of a non-conformal negative mode for the lower dimensional black hole 
solution, with some eigenvalue A < 0. This can be oxidized to give a negative mode of the black 
brane. This will satisfy the Euclidean versions of equations |4. 20| and [4.21| so it appears that 
this can now be Wick rotated and multiplied by exp(iniZ l ) (with fi 2 = —A) to obtain a classical 
Lorentzian threshold unstable mode. 



This conclusion is a little premature. Equations |4.20| and |4.21| are indeed recovered this way, 
but what about equations |4.22j and [4. 23| ? At first sight, it is hard to see how these could emerge 
from the negative mode since they are associated with the classical equations of motion in the 
flat extra dimensions, which the negative mode does not "know" about. However, if one takes 



the trace of equation 4.21 and combines it with equation 4.20, one can derive 



2V 2 Z + 



■ dZ + 2V fl Y tl 



2(3Y 



XZ. 



(5.12) 



The combination V u — f3d v cj) times equation |4.21| minus 2(k + f3 2 )d fJi (j) times equation |4.20| gives 
(after a long calculation involving the background equations of motion, the Bianchi identities 
and spherical symmetry) 



a, ( -v 2 z + poz ■ d<f> + v„y - (3Y ■ 



XY„. 



It follows that if A 7^ then 



Y, 



d^Z = 0, 



(5.13) 
(5.14) 
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i.e., 

V, (H^ - l -H?r v ) = PH^dJ + 2(k + f3 2 )fd^ + (3d J. (5.15) 

This looks very much like a gauge condition on H^. One is, of course, free to impose a 
gauge condition on since this is just the Euclidean metric perturbation for the black hole 
background, and no gauge condition has yet been imposed on this perturbation. However, it 
seems rather odd that a gauge condition should emerge from the equations of motion! The 
resolution of this puzzle is simply that, for A 7^ 0, the above equations are not equations of 
motion, and in fact the right hand side of the eigenvalue equation is not gauge invariant. A 
similar phenomenon would occur in Maxwell theory: the equation d v 8F^ v = \5A^ implies 
d ■ 5A = when A 7^ 0. Therefore, this eigenvalue equation is only valid Lorentz gauge. In a 
different gauge, the right hand side of the eigenvalue equation would be different. 



Further progress can be made by using equation |5.14| to eliminate from equation |5.12| , 
giving 

- V 2 Z + (3dZ ■ 00 = XZ, (5.16) 

which implies 

J d D x^e^\dZf = \J d D xy^e-^Z 2 . (5.17) 

The left hand side of this equation is positive. A normalizable negative mode would give a 
negative right hand side, so it follows that such a mode must have Z = 0. Using equation 5.14 
it follows that a normalizable negative mode must have Z = = 0, so equations [4.22| and 



4.23| are indeed satisfied. This completes the proof that that the threshold classical Lorentzian 



unstable mode in the gauge |4.19| is obtained from the negative mode in the gauge |5.15| simply 
by Wick rotation and multiplying by exp^/ijZ*) with fi 2 = —A. 



6 Conclusions 

It has been explained why there should be a relationship between classical and thermodynamic 
stability of black branes, as conjectured by Gubser and Mitra. Various assumptions were made 
which means that the above discussion falls well short of a rigorous proof. It is convenient to 
summarize the arguments used above, to emphasize where these assumptions were made. 

A) Local thermodynamic stability of a magnetically charged black brane is equivalent to 
positivity of the specific heat of the black hole obtained by dimensional reduction. If the black 
brane is locally thermodynamically unstable then so is the black hole, so it will have a negative 
mode. It was assumed that this is an s-wave, but this is probably acceptable since it seems 
very likely that higher partial waves would have higher eigenvalues. This negative mode can 
be oxidized to give a negative mode of the black brane. If a suitable gauge is chosen then this 
negative mode can be converted into a classical threshold unstable mode as described above. 
Such a mode would be expected to separate stable short wavelength fluctuations of the black 
brane from unstable long wavelength fluctuations, so the existence of such a mode indicates 
that a classical Gregory-Laflamme instability exists. So local thermodynamic instability implies 
classical instability. 

B) For the converse, if a black brane is classically unstable then the instability is expected 
to occur only at long wavelengths. Therefore there should exist a threshold unstable mode that 
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is independent of time. Provided this mode is of the longitudinal s-wave form discussed by 
Gregory and Laflamme, then it is possible to choose a gauge in which the components of the 
metric perturbation along the spatial worldvolume directions vanish. Dropping the exp(z/ijZ J ) 
factors and Wick rotating then gives an s-wave negative mode of the black hole. This implies 
the existence of a local thermodynamic instability of the black hole, assuming that the results 
of extend to the black holes considered here. It then follows that the black brane must also 
be locally thermodynamically unstable. 

The assumptions involved in proving A) that local thermodynamic instability implies clas- 
sical instability seem rather mild. The assumptions involved in proving the converse B) are less 
obvious, so it would be interesting to see if these could be given better justification. To see why 
proving B) rigorously is hard, consider the case when the black brane has positive specific heat. 
Proving B) in this case amounts to proving that the black brane is classically stable, which 
involves considering all fluctuations around the black brane, not just those of the GL form. 

If a black brane is sufficiently far from extremality then it should behave much like an 
uncharged black brane and therefore have negative specific heat. It then follows from argument 
A) above that the black brane will be classically unstable. As the brane approaches extremality, 
argument A) establishes that the brane will be classically unstable as long as the specific heat 
remains negative. So argument A) extends the results of GL to any magnetically charged black 
brane with negative specific heat. If argument B) is correct then the instability should vanish 
when the specific heat becomes positive. 

It has been explained that the reason GL found an instability that persisted all the way 
down to extremality is that the branes they studied always have negative specific heat. Had 
they chosen to work with branes whose specific heat became positive near extremality then, 
according to B), they would not have found an instability near extremality. 

For a black brane of a type II supergravity theory, with the exception of the D3-brane, it 
is always possible to dualize the field strength so that the brane is magnetically charged and 
therefore falls within the class of solutions considered in this paper. The non-extremal NS5, D5 
and D6 branes always have negative specific heat so A) implies that they are always classically 
unstable. However, the Fl, Dl, D2 and D4-branes have a specific heat that becomes positive 
near extremality. A) establishes that they should be classically unstable when their specific heat 
is negative, and B) implies that they should become classically stable near extremality. Note 
that the F-string and D-string have the same stability properties, as do the NS5 and D5-branes, 
as one would expect from the SL(2, R) symmetry of IIB supergravity. The major omission from 
this paper is, of course, the black D3-brane. I think that it is likely that the above analysis 
could be extended to cover this well as the other dyonic black brane solutions given in 

||26|| . If this is the case then the black D3-brane should be stable near extremality. 

The M2 and M5-branes are not covered by the analysis above but one can infer their stability 
properties from the behaviour of the F-string and D4-branes, which are obtained by dimensional 
reduction on a worldvolume direction. If the higher dimensional brane were classically unstable 
to a perturbation of the form considered above then one could reduce in a direction orthogonal 
to fii and the instability would persist in the lower dimensional solution. Since a GL instability 
is not expected for the near-extremal F-string or D4-brane, this implies that the near extremal 
M2 and M5-branes must be classically stable, in agreement with the specific heat analysis. 

The GM conjecture is supposed to cover any black brane with a non-compact translational 
symmetry. The branes considered in this paper were all assumed to have flat worldvolumes. 
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However, it might be possible to extend the analysis to branes with Ricci flat spatial worldvolume 
metrics of the "Kaluza-Klein" form 



ds 2 = e a *(dr + A) 2 + e^jdtfdyi, (6.1) 

where A is a 1-form independent of r, the coordinate along the non-compact direction. The 
classical instability would presumably involve fluctuations proportional to exp(z/xr). 

Perhaps the most striking result of this paper is that it is possible to show that many black 
branes are classically unstable without resorting to the arduous numerical analysis of |2j. 
Instead one just has to calculate the specific heat. So a simple thermodynamic calculation 
replaces a difficult classical calculation. This can be contrasted with black hole thermodynamics, 
where the temperature and entropy are, in principle, obtained by a long calculation [Q, but in 
practice are determined by simple classical quantities, namely the surface gravity and horizon 
area. Understanding these connections between classical and quantum properties of black holes 
remains one of the most interesting problems in quantum gravity. 
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